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IMPORTANT NOTICE 


Texas Instruments (TI) reserves the right to make changes in the 
devices or the device specifications identified in this publication 
without notice. Tl advises its customers to obtain the latest version 
of device specifications to verify, before placing orders, that the 
information being relied upon by the customer is current. 


Tl warrants performance of its semiconductor products, including SNJ 
and SMJ devices, to current specifications in accordance with TI's 
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testing of all parameters of each device is not necessarily performed. 


In the absence of written agreement to the contrary, Tl assumes no 
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under any patent right, copyright, or other intellectual property right 
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INTRODUCTION 


Algorithms, software, and hardware for designing and 
implementing a digital control system using the Texas 
Instruments TMS32010 signal-processing microprocessor are 
presented in this application report. Microprocessors, such 
as the TMS32010, are increasingly being used to implement 
algorithms for the control of feedback systems. The major 
factors contributing to this trend are increased availability 
and lower cost of suitable digital hardware. 

Current and potential applications include servo motor 
control, process control, robot arm and disk head controllers, 
and temperature and pressure controllers. Military and 
aerospace applications include stabilized platforms, flight 
control and autopilot systems, inertial reference systems, and 
general servomechanisms. 

In meeting control system requirements, designers face 
many alternatives. Cost, size, weight, power, and reliability 
decisions are typically application dependent. This report 
highlights the tradeoffs in developing algorithms, software, 
and hardware for a digital control system. 


DIGITAL CONTROL SYSTEMS 


General Considerations 

A digital controller is a signal processing system that 
executes algebraic algorithms inherent to the control of 
feedback systems (i.e., compensator and filter algorithms). 
Together with the plant (system to be controlled) and signal- 
acquisition circuitry, the digital controller makes up a digital 
control system such as the one shown in Figure 1. 

Note that the system requires analog-to-digital (A/D) 
converters for the external (command) inputs and for the 
state-variable feedback inputs to the digital controller. The 
system also requires a digital-to-analog (D/A) converter for 
the control outputs to the plant. 
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CONTROLLER 


The advantages of the digital control approach over the 
analog approach are: 

1. Ability to implement advanced control 
algorithms with software rather than special- 
purpose hardware 

2. Ability to change the design without changing 
the hardware 

3. Reduced size, weight, and power, along with low 
cost 

4. Greater reliability, maintainability, and testability 

5. Increased noise immunity. 


Microprocessor Selection and System 
Development Cycle 

Choosing an appropriate microprocessor is an 
important factor in efficiently implementing a digital control 
design. A class of special-purpose (as opposed to general- 
purpose) digital signal-processing microprocessors has been 
developed to enable fast execution of digital control 
algorithms. The Texas Instruments TMS32010 provides 
several beneficial features for implementing digital control 
system elements through its architecture, speed, and 
instruction set. 

A prominent feature of the TMS32010 is the on-chip, 
16x 16-bit multiplier that performs two’s-complement 
multiplication and produces a 32-bit product in a single 
200-ns instruction cycle. The TMS32010 instruction set 
includes special instructions necessary for fast 
implementation of sum-of-products computations 
encountered in digital filtering/compensation and Fourier 
transform calculations. Most of the instructions critical to 


- signal processing execute in one instruction cycle. References 


[1,2] give full details of the TMS32010 hardware and 
software considerations. 

Many system development tools are available and may 
be used for digital control system design.3 Figure 2 outlines 
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Figure 1. Digital Control System 


IDEA/SYSTEM SPECIFICATIONS 


DIGITAL CONTROL SYSTEM DESIGN PHASE, DISCRETE TRANSFER FUNCTIONS, 
STABILITY ANALYSIS, SIMULATION, ETC. 
(DIGITAL FILTER DESIGN AIDS FOR THE TMS32010) 


SYSTEM ARCHITECTURE DESIGN AND 
EVALUATION OF PROCESSOR PERFORMANCE 


(TMS32010 EVALUATION MODULE AND 
ANALOG INTERFACE BOARD) 













SOFTWARE CONFIGURATION 
AND DEBUGGING 
(TMS32010 ASSEMBLER/LINKER, 
SOFTWARE LIBRARIES, 
TMS32010 SIMULATOR) 


HARDWARE CONFIGURATION 
AND DEBUGGING 
(TMS32010 EMULATOR (XDS)) 


SYSTEM BINDING PHASE. HARDWARE/SOFTWARE INTEGRATION 
(TMS32010 EMULATOR (XDS)) 


Figure 2. Digital Control System Development Phases 


the system development cycle and ties the development tools 
to different project phases. 

Among the non-TI development aids, an interactive 
program called the Digital Filter Design Package developed 
by Atlanta Signal Processors Incorporated, may be useful 
in digital control system design. This program designs 
various types of filters, compensators, and other structures. 
It can be used when, for example, there is a need for several 
notch filters to filter out unwanted frequencies in a digital 
feedback loop. 


DIGITAL COMPENSATOR DESIGN 


Alternative methods exist for designing digital 
compensators. This section outlines several approaches to 
digital controller design and points out the analytical tools 
useful in the design process. 


Design Based on Analog Prototype 

A commonly used method of designing a digital control 
system is to first design an equivalent analog control system 
using one of the well-known design procedures. The resulting 
analog controller (analog prototype) is then transformed to 
a digital controller by the use of one the transformations 
described below. 

The design of the analog controller may be carried out 
in the s-plane using design methods such as root-locus 
techniques, Bode plots, the Routh-Hurwitz criterion, state- 
variable techniques, and other graphic or algebraic methods. 
The purpose is to devise a suitable analog compensator 
transfer function which is transformed to a digital transfer 
function. This digital transfer function is then inverse z- 
transformed to produce a difference equation that can be 
implemented as an algorithm to be executed on a digital 
computer. Two of the analog-to-digital transformation 


methods, the matched pole-zero and the bilinear 
transformation, are described as follows: 

1. The matched pole-zero (matched Z-transform) 
method maps all poles and zeroes of the 
compensator transfer function from the s-plane 
to the z-plane according to the relation: 


7. = esT 
where T is the sampling period. 


If more poles than zeroes exist, additional 
zeroes are added at z = —1, and the gain of the 
digital filter is adjusted to match the gain of the 
analog filter at some critical frequency (e.g., 
at DC for a lowpass filter). This method is 
somewhat heuristic and may or may not 
produce a suitable compensator. 

2. The bilinear (Tustin) transformation method 
approximates the s-domain transfer function 
with a z-domain transfer function by use of the 
substitution: 


z—1 
z+1 
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As in the matched pole-zero method, the 
bilinear transformation method requires 
substitution for s. Compensators in parallel or 
in cascade maintain their respective structures 
when transformed to their digital counterparts. 
This substitution maps low analog frequencies 
into approximately the same digital frequencies, 
but produces a highly nonlinear mapping for 
the high frequencies. 


To correct this distortion, a frequency 
prewarping scheme is used before the bilinear 
transformation. The frequency prewarping 
operation results in matching the single critical 
frequency between the analog domain and the 
digital domain. To achieve this result, the 
prewarping operation replaces each s in the 
analog transfer function with (wo/wp)s where 
Wo is the frequency to be matched in the digital 
transfer function and 


coe WoT 
Wp = mA tan —— 


Bilinear transformation with frequency prewarping 
provides a close approximation to the analog compensator? 
and is the most commonly used technique. Other methods 
for converting a transfer function from the analog to the 
digital domain are: the method of mapping differentials!8, 


the impulse-invariance method®, the step-invariance 
method!8, and the zero-order hold technique.® 

The basic disadvantage of design based on an analog 
prototype is that the discrete compensator is only an 
approximation to the analog prototype. This analog prototype 
is an upper bound on the effectiveness of the closed-loop 
response of the digital compensator. 


Direct Digital Design 

The direct digital design technique is a design method 
where a digital control system design is carried out from the 
very beginning in the z-domain to produce a digital control 
algorithm. Various pole-placement and relocation techniques 
are used to position the poles and zeroes of the system’s z- 
domain transfer function to yield the required system 
performance. | 

The root-locus method in the z-plane is similar to the 
root-locus design in the s-plane in that both are based on 
observing the position of the closed-loop poles as a function 
of the system gain. However, the effect of the locations of 
the poles and zeroes on system performance is not the same 
as in the s-plane. Knowledge of such correspondences in the 
z-plane allows those familiar with continuous system design 
to design digital compensators. For example, in stabilizing 
an unstable system, the adjustable gain is used to move the 
poles inside the unit circle instead of inside the left-half plane. 
A phase-lead controller may be employed to shift the root- 
locus to the left, which results in a system that responds 
faster. A phase-lag controller may be designed to allow for 
higher loop gain to produce smaller steady-state errors and 
improved disturbance rejection. !0 

The z-domain designer may also use pole-zero 
cancellation. In this technique, some of the poles and zeroes 
of the digital transfer function of the plant may be cancelled 
by zeroes and poles of the digital compensator. The 
compensator then introduces additional poles and zeroes at 
locations that enable the designer to achieve the desired 
performance characteristics. This method may affect system 
stability due to ‘‘inexact cancellation’’. If the poles of the 
plant that are to be cancelled lie close to the unit-circle, 
inexact cancellation may cause the system root-locus to go 
outside the unit circle at some point. This makes the system 
conditionally stable or unstable. 

In a system where fast response to the control input 
is required, a ‘‘deadbeat’’ approach may be taken. The 
deadbeat controller cancels all the zeroes and poles of the 
plant and introduces a pole at z = 1. The result is that the 
system output reaches its steady-state value in one sampling 
period with no overshoot. In practice, an ideal deadbeat 
controller is difficult to implement because of inexact 
cancellation. Although the output does not experience 
overshoot, it may oscillate between sampling instants. The 
design is ‘‘tuned’’ in the sense that the system response may 
be acceptable for a step input but not acceptable for other 
inputs. 

With z-plane design, conventional design techniques 
can be used to place the closed-loop system poles exactly 
where desired. The approximations associated with digitizing 


an analog prototype are thereby eliminated. The disadvantage 
of this technique is the relative difficulty of visualizing the 
effect of pole-zero locations in the z-plane on system 
performance. To overcome the disadvantage of the z-plane 
technique, the designer can use the w-plane or, better yet, 
the w’-plane design technique.> Both techniques transform 
the design to a plane similar to the s-plane by means of the 
same kind of substitution as in the bilinear transformation 
described earlier. This procedure thereby allows the use of 
the familiar s-plane and frequency-domain methods of 
continuous-system design. The designer proceeds by first 
transforming the continuous plant to the z-domain and thence 
to the w-plane or w’-plane. The appropriate compensator is 
then devised and transformed back to the z-plane, where it 
is used to specify the corresponding computational algorithm. 


State-Variable Design Methods 

State-variable design methods can also be used, 
including state-variable feedback and optimal control based 
on quadratic synthesis. 

In the state-variable feedback technique, all of the states 
are measured and fed back through constant gains. This 
allows all of the closed-loop poles to be positioned at any 
desired locations in the z-plane, but does not affect the 
positions of the system zeroes. 

In a design based on quadratic synthesis, a performance 
index or cost function is minimized by proper choices of the 
control law or feedback compensator. In the most practical 
design, the resulting compensator has the same form as that 
resulting from the application of direct digital design 
techniques.’ 


DIGITAL COMPENSATOR IMPLEMENTATION 


Digital compensator algorithms execute on processors 
that use finite-precision arithmetic. The signal-quantization 
errors associated with finite-precision computations and the 
methods for the handling of these errors are presented in this 
section. 


Fixed-Point Arithmetic and Scaling 

Computation with the TMS32010 is based on the fixed- 
point two’s-complement representation of numbers. Each 
16-bit number has a sign bit, i integer bits, and 15-1 fractional 
bits. For example, the decimal fraction +0.5 may be 
represented in binary as 


0.100 0000 0000 0000 


This is Q15 format since it has 15 fractional bits, one 
sign bit, and no integer bits. The decimal fraction + 0.5 may 
equivalently be represented in Q12 format as: 


0000.1000 0000 0000 
This number is in the Q12 format because it has 12 


fractional bits, one sign bit, and three integer bits. Note that 
the Q15 notation allows higher precision while the Q12 


3 


notation allows direct representation of larger numbers. 

For implementing signal-processing algorithms, the 
Q15 representation is advantageous because the basic 
operation is multiply-accumulate and the product of two 
fractions remains a fraction with no possible overflows during 
multiplication. When the Q12 format is used, a software 
check for overflow is necessary. The subsection, 
OVERFLOW AND UNDERFLOW HANDLING, provides 
a detailed analysis of overflow handling. 

In the case where two numbers in Q15 are multiplied, 
the resulting product has 30 fractional bits, two sign bits, 
and (as expected) no integer bits. To store this product as 
a 16-bit result in Q15, the product must be shifted left by 
one bit and the most-significant 16 bits stored. The 
TMS32010 instruction SACH allows for this one-bit shift. 

In the case where a Q15 number is to be multiplied 
by a 13-bit fractional signed constant represented as a Q12 
number, the result (to correspond with Q15) must be left- 
shifted four bits to maintain full precision. The TMS32010 
instruction SACH allows for the appropriate shift. The Q15 
and Q1i2 representations are used in the example in the 
section, DESIGN EXAMPLE: RATE-INTEGRATING 
GYRO STABILIZATION LOOP. 

When fixed-point representations are used, the control 
system designer must determine the largest magnitudes that 
can occur for all variables involved in the computations 
required by the digital compensator. (Floating-point 
representations allow larger magnitudes, but take more time 
for the microprocessor to perform the required 
computations.) Once these largest magnitudes are known, 
scaling constants can be used to attenuate the compensator 
input as much as necessary to ensure that all variables stay 
within the range that can be expressed in the given 
representation. 

Several methods are used to determine bounds on the 
magnitudes of the variables. One method, called upper-bound 
scaling, provides a useful, although sometimes too 
conservative, bound on the magnitude, yet it is 
Straightforward to calculate. Consider a variable y(n) that 
is obtained as the output of a digital compensator H(z) when 
the input is the sequence x(n). The bound on y(n) is given by 


oo 


| Xmax | ue 
n= 1 


|h(n) | 


Ymax = 


where Xmax is the maximum value in x(n) and the sequence 
h(n) is the unit-sample response sequence for the digital 
compensator H(z). 

Other methods for estimating the upper bound are Lp- 
norm scaling, unit-step scaling, and the averaging 
method.?,10 

After ymax is determined, the scale factor can be 
chosen as the multiplier that is applied to x(n) prior to the 
compensator computations to ensure that y(n) remains within 
the required bounds. In addition, the control system designer 
may have knowledge concerning the bounds on the 
compensator variables based on prior experience, the 


4 


characteristics of the corresponding variables of analog 
prototypes, and simulation results. 


Finite-Wordlength Effects 

All variables involved in the digital compensator — the 
input, the compensator coefficients, the intermediate 
variables, and the output — are represented as finite- 
wordlength numbers. This restriction gives rise to errors. 
Another source of errors is the truncation or rounding that 
takes place when the 32-bit product of two 16-bit numbers 
is stored as a 16-bit number. Both of these errors give rise 
to the finite-wordlength effects discussed in this section. 

The representation of the compensator input as a finite- 
precision (quantized) number produces an input-quantization 
error. The size of this error for a rounding scheme can be 
anywhere from —(2~8)/2 to (2—8)/2 where B is the 
number of bits in a word. The input-quantization error is 
usefully modeled as a zero-mean random variable uniformly 
distributed between its positive and negative bounds. A 
technique>.!0 is available to calculate the variance of the 
corresponding error at the compensator output (its mean is 
zero). In this manner, the designer can determine the effect 
of input quantization on the compensator output. 

Similar quantization errors are associated with the 
multiplication process. Each multiplication is assumed to 
produce the “‘true’’ product with an error that is a zero-mean, 
uniformly distributed random variable. The variance of the 
corresponding error at the compensator output can be 
calculated in the same manner as for the error due to input 
quantization. These individual variances are then added to 
measure the total effect at the compensator output for each 
truncation or rounding. 

Another way to describe the effects of truncation or 
rounding is in terms of ‘‘limit cycles’’ which are sustained 
oscillations in the closed-loop system. These oscillations are 
caused by nonlinearities within the loop. In this case, the 
nonlinear quantizations are associated with the 
multiplications. Limit cycles persist even when the system 
input goes to zero, and their amplitude can be sizeable. No 
general theory is available to treat this nonlinear 
phenomenon. Bit-level simulations which model the 
compensator and the complete closed-loop system are used 
to ascertain their presence and effect on the closed-loop 
performance. 

When a digital compensator is implemented as an 
algorithm to be executed on finite-precision hardware, a 
problem arises with implementing the coefficients present 
in the corresponding transfer function (see section, 
TMS32010 IMPLEMENTATION OF COMPENSATORS 
AND FILTERS). The infinite-precision compensator 
coefficients must be rounded and stored using a finite-length, 
fixed-point binary representation. Due to this coefficient- 
quantization effect, the performance of the implemented filter 
will deviate from the performance of the designed digital 
filter. | 

The deviation in performance can be estimated by 
computing the filter’s pole and zero locations and the 
corresponding frequency response magnitude and phase for 


the compensator with the quantized coefficients. Coefficient 
quantization forces the filter’s poles and zeroes into a finite 
number of possible locations in the z-plane and is of most 
concern for filters with stringent specifications, such as 
narrow transition regions. 

The designer must choose the filter structure least 
sensitive to inaccurate coefficient representation. The choice 
should be of a modular rather than a direct filter structure. 
For example, a higher-order filter should be implemented 
as a cascade or parallel combination of first-order and second- 
order blocks. The reason for this choice is the lesser 
sensitivity to coefficient variations of the roots of low-degree 
polynomials in comparison with high-degree polynomials. 
Several methods for selecting the filter structures least 
affected by coefficient quantization are available.? 

To quantitatively evaluate the effect of coefficient 
quantization on the position of the poles or zeroes of a digital 
transfer function, a ‘“‘root sensitivity function’’ can be 
computed .6 


Overflow and Underflow Handling 

Digital control system algorithms are usually 
implemented using two’s-complement, fixed-point 
arithmetic. This convention designates a certain number of 
integer and fractional bits. The fixed-point arithmetic 
computations may, at some point, produce a result that is 
too large to be represented in a chosen form of fixed-point 
notation (e.g., Q12). The resulting overflow, if untreated, 
may cause degraded performance such as limit cycles and 
large noise spikes at the filter’s output which may contribute 
to the system’s instability. The system must be able to recover 
from the overflow condition, i.e., return to its normal, 
nonoverflow state. 

Consider an example of the Q12 representation. The 
number 7.5 multiplied by itself gives the result of 56.25, an 
overflow in Q12. However, no hardware overflow occurs 
in the accumulator; 1.e., 


7.5 0111.1000 0000 0000 
x 7.5 _0111.1000 0000 0000 
56.25 0011 1000.0100 0000 0000 0000 0000 0000 


For the Q12 representation, the above 32-bit product 
is shifted left four bits and the left-most 16 bits are retained: 


1000.0100 0000 0000 


The correct answer is 56.25, but the number stored in 
the Q12 representation is — 7.75. 

The TMS32010 has a built-in overflow mode of 
operation that, if enabled, causes the accumulator to saturate 
upon detection of an overflow during addition when the 
accumulator register overflows. During multiplication, an 
overflow of the fixed-point notation may also occur even 
though the hardware overflow of the accumulator register 
does not occur. This is because the 32-bit result of a 
multiplication of two 16-bit numbers must be stored in a 


16-bit memory word in the form consistent with the chosen 
fixed-point notation (see above example). To adjust the 
location of the binary point, the storing operation requires 
that the number in the accumulator be shifted left and 
truncated on the right before storing. If the most significant 
bits shifted out contain magnitude information in addition 
to sign information, an overflow in the chosen fixed-point 
notation results. 

To track overflows associated with the number 
representation, the control system software should contain 
an appropriate overflow-checking routine in those places 
where multiplications and additions occur. This routine 
should not rely exclusively on the TMS32010’s overflow 
mode to intercept and correct the overflow occurrences. 

Two approaches may be used to handle overflows. The 
first is to prevent the overflow from occurring by choosing 
conservative scaling factors for the numbers used in 
computations, as described in the subsection, FIXED-POINT 
ARITHMETIC AND SCALING. These scaling factors are 
used to limit the range of inputs to each of the basic building 
blocks of the compensator, namely, the first- and second- 
order filter sections. The scaling factor chosen reduces the 
input magnitude and consequently all other signal levels, 
thereby enabling the compensator coefficients, the expected 
inputs, and their products and sums, all to be represented 
without overflow. The scaling must also maintain the signal 
levels well above the quantization noise. 

The second approach for handling overflow is to adjust 
the sum or product each time an overflow occurs. To 
accomplish this, an overflow checking routine must be 
written and executed at certain points along the computational 
path. The routine must check whether the number just 
computed and residing in the 32-bit accumulator can be stored 
without overflow in a 16-bit memory location in accord with 
the chosen fixed-point notation. Once the routine detects an 
overflow condition, it should replace the computed number 
with the maximum or minimum representable two’s- 
complement number. This scheme simulates a saturation 
condition present in analog control systems. To prevent 
overflow limit cycles, the saturation overflow characteristic 
is preferred to the two’s-complement, ‘‘wrap-around’”’ 
characteristic .? 

An example of the overflow checking and correcting 
technique for a first- and second-order filter subroutine is 
provided in the section, TMS32010 IMPLEMENTATION 
OF COMPENSATORS AND FILTERS. This Direct-Form 
II implementation subroutine checks for overflow 
occurrences upon computation of the filter’s intermediate 
state variable and again upon computation of the filter’s 
output. 

In a digital control system, the first- and second-order 
building blocks are either cascaded or connected in parallel 
to compute a series of control algorithms. The first- and 
second-order filter subroutine, called to compute each of the 
control system’s elements, uses 16-bit memory locations as 
storage media for its intermediate values, in which case it 
is appropriate to check for overflow in each block. 


At the end of a computational chain — before the final, 
computed digital output is ready for transfer to the analog 
domain — it is necessary to check that the number being sent 
to the digital-to-analog converter is within the range based 
on the manner in which the converter is interfaced to the 
processor data bus. For example, if a 12-bit converter is 
wired to the 12 least significant bits (LSBs) of the 16-bit 
processor data bus, then the 12 LSBs must contain both 
magnitude and sign information, which may require that the 
original 16-bit number be adjusted or limited before being 
sent to the converter. 

Underflow conditions, which can also appear during 
digital control algorithms, are conceptually similar: to 
overflows in that the computed value contained in the 32-bit 
accumulator is too small to be accurately represented in a 
16-bit memory word in the chosen fixed-point notation. One 
possible solution to this problem is to multiply the small result 
by a gain constant to raise its value to a representable level. 
The appropriately chosen gain constant may come as a result 
of gain distribution throughout the digital control system, 
whereby large gains from some of the building blocks get 
uniformly distributed over a range of the system’s sections. 


DIGITAL CONTROL SYSTEM 
SOFTWARE DESIGN PHILOSOPHY 


To maximize the manageability and portability of the 
system software, a modular or top-down design technique 
should be used. This section shows how the modular software 
structure and the proper layout of system memory contribute 
to the efficient implementation of a digital control design. 


Modular Software Structure 

The concept of modular software design is a technique 
developed to make system software more manageable and 
portable. Top-down design is used to break up a large task 
into a series of smaller tasks or building blocks, which in 
turn are used for structuring a total system in a level-by-level 
form. At the end of a top-down design process, a number 
of modules are linked together which, under the control of 
a main program, perform as a complete system. 

In addition to making the software-development and 
software-modification processes more manageable, modular 
design also enhances software portability. Digital control 
systems use a number of standard functional blocks such as 
compensators, notch filters, and demodulators. It is therefore 
likely that a designer who already has access to one digital 
control system will want to ‘‘borrow’’ some of its functional 
building blocks to quickly implement a new, different control 
unit or reconfigure the existing one. The designer who has 
access to these functional blocks or modules needs only 
modify the main program by providing a different sequence 
of subroutine calls. An initialization routine, a first- and 
second-order filter routine, a roundoff routine, and an 
overflow checking routine are examples of functional 
building blocks. 


Each software module is written as a subroutine with 
a clear and efficient interface (for parameter passing, stack 
use, etc.) with the main program. In order to maintain the 
general-purpose function of the module, the data used in 
computations within a module (i.e., filter coefficients, state- 
variable values, etc.) should be accessed using indirect 
addressing rather than direct addressing. Only those variables 
whose values remain unchanged should be addressed directly. 


Layout of TMS32010 Data Memory 

The layout of the TMS32010 data memory in a digital 
control system implementation should be defined in 
accordance with the requirements of the software modules 
used in the implementation of the system. The procedure is 
illustrated by the first- and second-order filter subroutine of 
the section, TMS32010 IMPLEMENTATION OF 
COMPENSATORS AND FILTERS. This subroutine 
manipulates its pointer registers so that upon completion of 
the computations in one filter section, the registers 
automatically point to the set of coefficients and state 
variables of the next filter section. 

If the software designer arranges his filter coefficient 
and state-variable sections in the order of execution of the 
control-system algorithms, a sequence of compensators and 
filters may be executed with a single subroutine call for each 
element. This scheme enables faster execution of the control 
algorithms since there is no need to explicitly reload the 
pointers in order to match the requirements of the current 
software module being called. 

A designer must define all of the data memory locations 
(set up a system memory map) at the beginning of the 
program. An efficient way to accomplish this is to use the 
TMS32010 assembler’s DORG (dummy origin) directive. 
This directive does not cause code generation. DORG defines 
a data structure to be used by the system; 1.e., it generates 
values corresponding to the labels of consecutive data 
memory locations. Using the DORG directive, as opposed 
to equating labels with data memory locations through the 
EQU directive, provides flexibility when the data structure 
needs to be modified. For example, when defining a number 
of new data memory locations, the labels are inserted in the 
middle of the ‘‘dummy’’ block and the assembler assigns 
the values automatically. This function would have to be 
performed manually if the EQU directive were used. 

The software designer must also build a table in the 
TMS32010 program memory that corresponds to the 
previously defined data memory map. The table is then 
loaded into the data memory by the initialization routine 
during system startup. 

These techniques are illustrated in the next section, 
TMS32010 IMPLEMENTATION OF COMPENSATORS 
AND FILTERS. Note that in the example program in 
Appendix B, location ONE has to be the last location in the 
table. Note also that the states and the coefficients of the 
filters are defined in reverse order to the order in which the 
filters execute. This is due to the way the initialization and 
filter routines are written. 


TMS32010 IMPLEMENTATION OF 
COMPENSATORS AND FILTERS 


Design procedure and error handling for the standard 
first- and second-order compensator and filter subroutine are 
described in this section. Methods for implementing higher- 
order structures, implementation tradeoffs, and examples of 
typical compensators and filters are also given. 


Standard First- and Second-Order Block 
as a Subroutine 

A standard first- and second-order compensator section 
is a prime example of the building block philosophy discussed 
earlier. The routine presented here computes first- and 
second-order IIR filter sections using the Direct-Form II 
network structure!2 and performs roundoffs and overflow 
checking. The Direct-Form II, although it somewhat obscures 
the definition of the variables, is chosen over the Direct-Form 
I because it requires fewer ‘“‘delays’’, 1.e., data storage 
locations, in its computational algorithm. 

Consider the second-order transfer function: 
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For Direct-Form II, the corresponding difference 
equations are: 


x(n) = e(n) — D1 xMm-—1) — D2 x(n-2) 


u(n) = NO x(n) + N1 x(n—1) + N2 x(n—2) 
The signal flowgraph for this transfer function is shown 
in Figure 3. 
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Figure 3. Direct-Form II Compensator/Filter 


The filter routine accommodates a scaling scheme as 
defined on the main program level; i.e., the values can be 
scaled by 215 or 212. The routine is written so that chain 
implementation of a number of compensators and filters is 
possible if the data structure, namely, the coefficient and 
State-variable tables, are properly arranged (see the previous 
section, CONTROL SYSTEM SOFTWARE DESIGN 
PHILOSOPHY). The routine takes its input from the 
accumulator and outputs the result to the accumulator so that 
the main program can efficiently call for the successive 
execution of the filter routine with a different set of 
parameters each time. 


Two checks for overflow are made within the filter 
routine. One is made upon computing the value of the 
intermediate (state) variable and the other upon computing 
the filter’s output. Each overflow check determines whether 
the 32-bit computed result can be stored in a 16-bit memory 
location under the adopted scaling scheme. If an overflow 
condition occurs, the routine “‘saturates’’ the output; 1.e., 
it returns the maximum or minimum representable value. 

A drawback of overflow checking upon computing the 
output is the loss of precision of the least-significant bits of 
the accumulator, which are truncated during the accumulator 
storing operation. This loss of precision is insignificant, 
however, in comparison with the loss of precision due to an 
overflow condition. 

The first-order filter section is computed exactly like 
the second-order section. The two coefficients N2 and — D2 
that multiply the ‘‘oldest’’ value of the intermediate (state) 
variable, i.e., the bottom branch of the filter in Figure 3, 
are equal to zero. This scheme reduces the second-order 
digital filter to the first-order filter. An example program 
that uses the first- and second-order filter routine to compute 
several elements of a digital control system is given in 
Appendix B. 


Higher-Order Filters: Cascade Versus 
Parallel Tradeoffs 

A higher-order filter or compensator in a digital control 
system can be implemented either as a single section or as 
a combination of first- and second-order sections. The single 
section or direct implementation form is easier to implement 
and executes faster, but it generates a larger numerical error. 
The larger error occurs because the long filter computation 
process involves a substantial accumulation of errors resulting 
from multiplications by quantized coefficients and because 
the roots of high-order polynomials are increasingly sensitive 
to changes in their (quantized) coefficients. For this reason, 
the direct realization form is not recommended except for 
a very low-order controller. 

The suggested method of implementing a high-order 
transfer function is to decompose it into first-order blocks 
(to accommodate single real poles) and second-order blocks 
(to accommodate complex conjugate poles or pairs of real 
poles), and connect these blocks either in a cascade or a 
parallel configuration. 

For the cascade realization (see Figure 4), the transfer 
function must be decomposed into a product of first-order 
and second-order functions of the form: 


D(z) = K Dj(z) Daz). . . Dp(Z) 
Each second-order block has the form: 


Dig) = NO + Ni z~! + N2z~-2 
mm 1 + Di z7-l + D227-2 
Each first-order block is obtained by equating the 
coefficients of (z~2), i.e., N2 and D2, to zero. 
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Figure 4. Cascade Implementation of a High-Order 
Transfer Function 


The designer must decide how to pair the poles and 
zeroes in forming the D;(z). A pole-zero pairing algorithm 
that minimizes the output noise is available.!9 The ordering 
of the D,(z) also affects output noise due to quantization and 
whether or not limit cycles are present. 10 

For the parallel realization (see Figure 5), the transfer 
function is expanded as a sum of first-order and second-order 
expressions of the form: 


D(z) =K + D,(z) + Do(z) +... + Dplz) 


where the first-order blocks have the form: 


____ NO 
Pe) Ee 
and the second-order blocks have the form: 
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Figure 5. Parallel Implementation of a High-Order 
Transfer Function 


The concept of ordering of Dj(z) does not apply to the 
parallel configuration. 

Pole-zero pairing is fixed by the constraints imposed 
by the partial-fraction expansion. 

The parallel realization has an obvious advantage over 
the cascade form when an algorithm executes on a 
multiprocessor system where the filter algorithm can be split 
up among the processors and run concurrently. 

No significant difference is apparent between the 
parallel and cascade realizations in performance factors such 
as execution speed and program/data memory use when the 
algorithms execute on a single processor.!3 The parallel 
algorithm could possibly provide more precision in 
computing the filter’s output if the designer decided to save 
the double-precision (32-bit) results from each of the first- 
and second-order sections and perform a double-precision 
addition to calculate the final filter output. 


Examples of Typical Compensators 
and Filters | 

Examples of structures used in digital control systems: 
compensators (lowpass filters), and notch filters are shown 
in Table 1. The analog and digital versions of the transfer 
functions, along with the scaled form of the transfer functions 
and their coefficient values, are given. A complete example 
of a transformation from the analog to the digital domain 
using bilinear transformation with frequency prewarping is 
presented in Appendix A. 

The gain constant that appears in some transfer 
functions can be implemented either by integrating it into 
its transfer function or by distributing it over a number of 
filter sections in a cascade implementation scheme. 


PROCESSOR INTERFACE CONSIDERATIONS 


Alternatives should be considered when designing the 
data-acquisition portion of the digital controller hardware. 
This section addresses the A/D and D/A converter selection, 
different analog sensor interface methods, and 
communication with the host processor. 


Table 1. Examples of Analog and Digital Versions of Common Transfer Functions 


Digital Transfer Function 


Analog Prototype Transfer Function (Computed by Bilinear 


Transformation with Frequency 
Prewarping; f, = 4020 Hz} 


First-Order Compensator: 


= 
D(z) = 0.012436 fee 


G(s) = 100 Se 
S+1 1.0-0.99975 z—1 


Second-Order Compensator: 


S2 +68.25 +3943 


G(s) = 1000 = 
ie) $2 +2512$+6.31 x 108 


D(z) = 870.77 


100-Hz Notch Filter: 


aie S2 + 3.9478 x 105 . desis 
$2 + 125.6648 + 3.9478 x 105 = 


1.0-1.9824 z-14+0.9826 z-2 
1.0-1.2548 z- 140.5474 z~2 


0.98467 — 1.94534 z—1+0.98467 z~2 
1.0-1.94534 z~1+0.96935 z~2 





Scaled Coefficients and 
Order of Storage in 
Data Memory (as required by 
First- and Second-Order 
Filter Routine) 

NO = 408 
N1 = 408 
D1 = 32760 


Scaled Digital Transfer Function 


Scaling Factor = 215 

408 + 408 z~! 
D(z) = as are ae Ir 
32768 — 32760 z-1 


Scaling Factor = 212 


4096 —8120 z~1+4025 z~2 


D(z) = 870.77 
” 4096-5140 z-1+2242 2-2 


4033 
~ 7968 
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7968 
D2 = —~3970 


Scaling Factor = 212 


4033 — 7968 z~1+4033 2-2 
4096 - 7968 z~1+3970 z-2 
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A/D and D/A Conversions and 
Integrated Circuits 

The A/D and D/A converter selection for a control 
system design may be based on several factors. Among the 
most crucial factors are the maximum conversion speed of 
the converter and the wordlength of the device. 

The A/D conversion speed relates directly to the 
required sampling rate of the specific application. This rate 
is determined by the need to sample fast enough to prevent 
aliasing and excessive phase lag and to sample slow enough 
to avoid the unnecessary expense and accuracy of high data 
rates. 

The A/D wordlength should be chosen based on a 
worst-case analysis using the following two criteria: 

1. The dynamic range of the continuous input 
signal, and 
2. The quantization noise of the A/D converter. 

For dynamic range, the designer should determine the 
minimum and maximum values of the continuous input that 
need to be accurately represented and select the A/D 
wordlength in bits based on the resolution required within 
this range. 

Quantization noise is due to the quantization effect of 
the A/D. The value of this noise during a single conversion 
can be represented by the difference between the exact analog 
value and the value allowable with the finite resolution of 
the A/D. This quantization noise may assume any value in 
the range —q/2 to +q/2 for a rounding converter or 0 to q 
for a truncating A/D converter where q is the quantization 
level. The quantization level q is equal to the full-scale 
voltage range divided by 2B where B is the number of bits 
in the converter. The quantization noise may be modeled as 
uniformly distributed noise. The designer should make his 


choice based on the maximum acceptable quantization level. 

The D/A converter wordlength should be chosen in a 
similar manner to choosing the A/D wordlength by 
considering the dynamic range of the output signal. 

The effects of A/D and D/A converter wordlength on 
the performance of a high-speed control system are detailed | 
in the University of Arkansas study (see the section on the 
design example of the TMS32010-based rate-integrating gyro 
positioning system). The study analyzed the time-domain 
performance of the system (unit impulse, step, ramp, and 
torque-disturbance response) as a function of A/D and D/A 
wordlengths. Twelve-, fourteen-, and sixteen-bit converters 
were used. The only significant difference found between 
them was the steady-state error. Twelve-bit converters were 
found to be adequate. 

In a multi-input digital control system, the signal 
acquisition portion of the digital controller must provide for 
the multiplexing of several analog inputs into a single A/D 
converter. Consequently, some external devices are needed 
to pre-filter (antialiasing filters), sample and hold the analog 
signals from each channel (S/H circuits), and multiplex the 
signals onto the A/D converter (analog multiplexer). 
Multiplexing and filtering may also be necessary at the output 
in cases where the digital control system computes multiple 
outputs for the control of the plant. 

Two configurations of a cost-effective, multichannel 
data acquisition system for a digital controller are shown in 
Figure 6. The first accommodates up to eight inputs; the 
second can accommodate up to 32 inputs. Note that in these 
two systems, only one S/H per eight inputs exists, and the 
variables are Sampled in sequence with the same sampling 
interval between successive samples of a given signal. There 
will be a ‘‘skew’’ in time between the samples of the various 
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Figure 6. Cost-Effective Data Acquisition Systems 


inputs with a possible unwanted effect on the system 
performance. In this case, the use of a fast A/D converter 
may be justified to minimize this effect. 

If truly simultaneous sampling is required, an array of 
S/H circuits may be used to capture the values of all the inputs 
concurrently. This solution, shown in Figure 7, is more 
expensive due to the cost of S/Hs. In such simultaneous 
sampling systems, a fast conversion must be performed 
before the signal values present on the S/Hs start to droop. 
Therefore, the maximum conversion rate must be fast enough 
to accommodate this constraint. 
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Figure 7. Data Acquisition System with 
Simultaneous Sampling 


In some cases, high-speed A/D converters (100- to 
500-kHz conversion rates) are required. Two 12-bit A/D 
devices that will accommodate these speed requirements are 
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the ADC85 and the AD5240 from Analog Devices. The 
ADC$835 allows a conversion rate of up to 100 kHz and the 
AD5240 up to 200 kHz. There are other converters available 
from several manufacturers. | 

D/A converters, on the other hand, are inherently 
faster, the selection is much broader, and the cost is less. 

When high-speed converters are not needed (when only 
a few input channels exist or a single converter per channel 
is justified), devices such as the TCM2913 and TCM2914 
codecs may be useful in digital control applications. Although 
telecommunications-oriented, these devices are low in cost 
and provide on-chip antialiasing and smoothing filters. The 
TCM2913 or TCM2914 both contain A/D and D/A 
converters. The 8-bit digital output of the A/D and the 8-bit 
digital input to the D/A are both arranged in a companded 
(compressed/expanded) form using yp-law or A-law 
companding techniques. The p-law and A-law companding 
techniques allow small numbers to be represented with 
maximum accuracy, but require a conversion routine before 
the companded samples can be used in two’s-complement 
computations. Such conversion routines are based on lookup 
tables and need only a few TMS32010 instruction cycles to 
execute.!4 The devices interface to the processor in a serial 
form and convert the data at a maximum rate of 8 kHz. 

All of these data acquisition systems can accommodate 
differential inputs from analog transducers, such as pressure 
sensors, strain gauges, and others. To maintain accuracy in 
the case of a low-level input signal and to minimize noise 
effects, twisted-pair leads can be used to connect the 
transducer output to an instrumentation amplifier 
(differential-to-single-ended conversion circuit) that in turn 
is connected to the analog multiplexer. Alternatively, 
balanced twisted-pair leads can be connected to a differential 
analog multiplexer which drives an instrumentation amplifier 
of the same kind. The amplifier rejects the common-mode 
noise and presents the single-ended output to the S/H circuit 
and the A/D converter.!5,16 These two configurations are 
shown in Figure 8. 
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Figure 8. Differential Input Configurations 
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Synchronization of the Processor 
and External Devices 

In a multichannel data acquisition system with one A/D 
converter, a designer must generate a sequence of timing 
signals to synchronize the S/H circuits, the analog 
multiplexer, the A/D converter, and the input/output latches 
with the operation of the TMS32010. The designer may 
decide to generate the timing signals from the TMS32010 
clock by subdividing its frequency or using a timing and 
control circuit based on its own clock. 

_ An alternative way to build a multi-channel data 
acquisition system is to designate a separate A/D for each 
channel. In this case, the timing-signal generation is simpler 
and the A/D converters used may be slower and less costly, 
although more of them are necessary. The designer should 
perform a tradeoff analysis based on board space, overall 
system cost, and power consumption. 


Communication with Host Computer 

In addition to having a fast signal-processing 
microprocessor, a need may exist for an executive processor 
to monitor the system’s operation. Such an executive 
processor would be used for system startup/initialization 
(coefficient and initial-condition loading), responding to 
emergency conditions such as overflow and underflow, 
system reprogramming/reconfiguration (loading a new 
program or a new set of coefficients), and a thorough system 
test and calibration. The system should be constructed so that 
the executive processor can interrupt, halt, or alter the 
execution of the signal processor at any time in response to 
contingency situations. 


DESIGN EXAMPLE: RATE-INTEGRATING 
GYRO STABILIZATION LOOP 


An example of using the TMS32010 processor to 
implement a digital control system is presented in this section. 
Sampling rate selection, the system’s hardware and software, 
and system performance are discussed. 


System Description 

The system used as an example of the application of 
the TMS32010 is a servo-control system for stabilizing a 
large, two-axis gimbaled platform with a DC-motor drive. 
Inertial rate-integrating gyroscopes mounted directly on the 
platform serve as angular motion sensors. Such systems are 
required for the precise control of line-of-sight (LOS) and 
line-of-sight rate for use in pointing and tracking applications 
for laser, video, inertial navigation, and radar systems. 

At present, digital control is not normally used in 
systems of this type because of the fast throughput rates and 
computational accuracy required to perform the control 
computations and notch filtering. Current line-of-sight 
stabilization systems continue to use analog electronics to 
implement servo-compensation functions and error-signal 
conditioning. Thus, the system is representative in 


complexity and performance of typical systems currently in 
use by the aerospace industry and are candidates for 
microprocessor-based digital control. 

The digital control system was designed as part of a 
research contract carried out by the University of Arkansas 
under the sponsorship of Texas Instruments from February 
1982 to February 1984. 18,19 


System Model and Control Compensation 

A single axis of the stabilization system has two primary 
control loops: the rate loop and the position loop. In addition, 
a tachometer loop exists within the position loop. The rate 
and position loops are identified in Figure 9, a diagram of 
the elevation axis of the system. In its analog version, the 
system employs analog electronics to implement all control 
compensation and signal conditioning functions. 

Figure 10 identifies those filters and compensators in 
the rate loop that are to be incorporated into the digital control 
system. 

This study’s approach provides a digital implementation 
of the designated analog elements of the rate loop without 
sacrificing closed-loop performance. In keeping with the 
recommendations of the DIGITAL COMPENSATOR 
DESIGN section, the technique for the conversions of the 
analog compensators and notch filters to their digital 
counterparts is the bilinear transformation with frequency 
prewarping. Within the rate loop, the transfer functions to 
be implemented digitally consist of a first-order and a second- 
order compensator, along with six notch filters. Within the 
position loop, there is one first-order compensator and one 
notch filter. The transfer functions, shown in Table 2, list 
both the analog prototypes and their digital equivalents. 

The sampling rate chosen is 4020 samples per second 
(sampling period is 249 ys). This rate is more than twice 
the highest frequency of consequence (1800 Hz, the highest 
rate-loop notch frequency) to prevent aliasing. The rate is 
fast enough to prevent excessive phase lag in the rate loop 
and is more than ten times the closed rate-loop bandwidth 
(approximately 80 Hz).° The rate was also chosen to be an 
integer multiple of 30 Hz, which is a commonly used update 
rate of the video and infrared imaging/tracker devices that 
provide the line-of-sight rate command to the stabilization 
system’s rate loop. The update rate of the imaging device 
and the sampling rate within the rate loop are thus 
synchronized. 

After simulating the closed rate loop, the phase margin 
was found to be five degrees less than it was for the all-analog 
system, due to the computational and other delays associated 
with sampling. To overcome this deterioration in phase 
margin, the second-order rate-loop compensator was 
redesigned to provide additional phase lead. The compensator 
was modified to provide enough additional phase lead so that 
the phase margin of the digital system matched that of the 
analog system. The modified compensator is listed in the 
table. 
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Figure 10. Line-of-Sight Stabilization/Pointing System Rate-Loop [19] 


Table 2. Analog and Digital Compensators and Notch Filters 


Compensator/ Auslow Teanater Eanciion Digital Transfer Function 
Filter Element ? (fs = 4020 Hz) 

























Rate-L 

oe G6 = 100 D6 = 0.1244 (1.0+1.0 Z~-1) 
- S+1 : 1.0-0.99975 zZ-1 

Compensator 











Rate-Loop 
2nd-Order 
Compensator 











754.7101 (1.0-1.98426 Z2-1+0.9845 2-2) 
1.0-1.255 Z-1+0.5474 2-2 


1000 (S2+62.8S + 3943) 
S2 + 2512S +6.31x106 












S2+ (22 1800)2 


G8 = 27 1800 
S24 ——— S+(2m 1800)2 










Rate-Loop 
1800-Hz 
Notch Filter 











0.96877 + 1.83411 Z-1+0.96877 Z-2 


D8 = 
1.0+1.83411 Z-1+0.93754 2-2 










S24 (2x2 900)2 


G9 = 2m 900 
S24 aoe S +(22 900)2 






Rate-Loop 
900-Hz 
Notch Filter 







0.8352 —-0.27291 Z~1+0.8352 Z-2 


D9 = 
1.0-0.27291 Z-1+0.67041 Z-2 














S2+ (22 560)2 


G10 = 27 560 
Si S +(2r 560)2 






Rate-Loop 
560-Hz 
Notch Filter 







0.9287 —- 1.19021 Z-1+0.9287 Z~-2 


D10 = 
1.0-1.19021 Z-1+0.8574 2-2 
















Rate-Loop $24 (2% 140)2 


140-Hz Gil = © 2n 140 
S24 S+(2r 140)2 
Notch Filter 5 







0.97875 — 1.91083 Z-1+0.97875 2-2 


D11 = 
1.0- 1.91083 Z-1+0.95751 Z-2 








S2+(22 120)2 


G12 = 22 120 
$24 S+(2x 120)2 


Rate-Loop 
120-Hz 
Notch Filter 










0.9817 — 1.92896 Z-1+0.9817 2-2 
1.0-1.92896 Z~1+0.96339 Z-2 







D12 = 











S2+(22 100)2 


G13 = 22 100 
S24 a S+(2x 100)2 






Rate-Loop 
100-Hz 
Notch Filter 









0.98467 — 1.94534 Z- 14+0.98467 Z-2 


D13 = 
1.0-1.94534 Z-14+0.96935 2-2 















Position-Loop 
1st-Order 
Compensator 










6.6S + 45.54 
Ss 


6.60566 —- 6.59434 z-1 
1.0-Z~-1 






D22 = 














S24 (22 400)2 
G24 = 2x 400 
5 






Position-Loop 
400-Hz 
Notch Filter 







0.94471 — 1.53204 Z-1+0.94471 2-2 
1.0-1.53204 Z-1+0.88942 2-2 








S+(22 400)2 
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Hardware SO - | 

In the digital control system, the analog compensators 
and the notch filters are replaced by a digital signal processor, 
the TMS32010, along with the additional interface hardware 
needed to provide the digital input signals to the controller 
and the analog signals to the plant. Figure 11 shows the 
system hardware block diagram. :. 

The hardware was packaged onto five wirewrap 
boards. It was fabricated as a prototype test bed, and was 
constructed from commercially available components that 
have military-specification counterparts. Twelve-bit A/D and 
D/A converters were used, based on the studies of time- 
domain performance characteristics of the system. 


Software | 

The TMS32010 software 1s composed of four modules: 
Initialization Routine, Main Program, Rate-Loop 
Subprogram, and Subroutines. Figure 12 shows the system 
software block diagram. 

The Initialization software disables and enables 
interrupts, loads data memory with filter coefficients, 
program constants, and gain terms, and initializes the 
TMS32010 registers. 

The Main Program software calls the Delay Subroutine 
at the beginning of each sample period to wait for the A/D 
to complete conversion of all input variables. It then does 
on-line compensation for the error signal sensor variation 
by executing the A/D Drift Subroutine. The Main Program 
then reads the value of the input variable, calls the Rate-Loop 
Subprogram to compute the control output, and, when that 
subprogram returns the output variable, loads it into the 
appropriate output register. | 

The Rate Loop Subprogram calls subroutines that 
perform each compensator and notch filter computation and 
checks the computed output for overflow. | 

The Subroutines consist of a single routine for 
performing any of the compensator or notch-filter 






GUARD 
FILTER 







SAMPLE 
& HOLD 
AMPLIFIER 
ARRAY 







8-CHANNEL 
MULTI- 
PLEXER 


ANALOG 
INPUTS 


A/D 
CONVERTER 






CONVERSION 
COMPLETE 





E INPUT 
4020-Hz SAMPLE LATCH 


AID START CONVERSION eon 


computations (first- and second-order filter routine), along 
with routines for checking overflow, providing delay, and 
performing multiplication of low-precision numbers by a 
constant. a 

The A/D Drift Subroutine compensates on-line for the 
variations in the rate-loop error signal sensor (as a function 
of time and temperature). The subroutine uses an external 


calibration input and follows the model of the sensor 


variations to estimate the true value of the A/D input. 

The digital control system is interrupt-driven. An 
interrupt occurs every 1/4020 seconds (approximately 250 
us). This starts the A/D conversion of a new set of sample 
inputs and restarts the TMS32010 on a new pass through its 
software. 

A TMS32010 Evaluation Module (EVM) and Emulator 
(XDS) were used in the software development to permit 
single-step execution of the software for comparison with 
the corresponding computations produced by simulations 


written with the aid of the Continuous System modeling 


Program (CSMP). These simulations take into account the 
input/output signal quantization levels, microprocessor 
architecture, memory and internal register lengths. 
Other software functions associated with a complete, 
self-contained control module include: 
1. System calibration, testing, and startup 
2. Error checking and contingency responses 
3. Setting of gains, time constants, and other 
programmable or adjustable parameters 
4. System shutdown. 
These functions are implemented by a general-purpose 
executive processor (SBP9989), thus allowing the TMS32010 
to handle computation-intensive tasks. 


System Performance 
The system performance was evaluated in the following 
two-step procedure: | 
1. A hybrid computer system was constructed 


TO SBP9989 BUS 


LI 


BUS TRANSMITTERS/ 
| RECEIVERS 
;~— i ct OUTPUT 
S. D/A ANALOG 
rmss2o10 ARRAY CONVERTER)" OUT 










OUTPUT 
LATCH 
CONTROL 


Figure 11. Digital Controller Hardware Block Diagram [19] 
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Figure 12. Digital Controller Software Block Diagram [19] 


consisting of an analog-computer implemen- 
tation of part of the nondigital portion of the 
rate loop coupled with the TMS32010-based 
digital controller. 

A full-scale CSMP simulation of the entire rate 
loop was conducted. 

The closed-loop performance of the rate loop was 
characterized by the following responses: rate-command step 
response, torque-disturbance step response, and torque- 
disturbance frequency response. The results are shown in 
Figures 13 through 15. 

In the rate-command step response, the percent 
overshoot, peak time, and settling time are similar for both 
the discrete and continuous systems, but the continuous 
system is slightly smoother. The torque-disturbance step 
response shows that the discrete system is slightly slower 
in correcting for a torque disturbance input. In addition, the 
discrete system has a low-level oscillation (limit cycle). The 
frequency responses for a torque-disturbance input are also 
similar, with the continuous system having slightly better 
torque disturbance rejection in the low-frequency region. 
These results show that the analog and the digital systems 
are comparable even though no special efforts were made 


2 


to take advantage of the capability that digital control offers. — 


The flexibility of the digital control system was 
demonstrated by programming the digital system with the 
capability to correct for a variation in the sensor input to the 
A/D converter. The system was able to correct on-line (by 
using a known standard, calculating the gain, and dividing 
it out) for a 50 percent sinusoidal variation in the sensor gain. 

The conversion between two different stabilization 
systems serves as another flexibility example. The software 


of a small, two-axis stabilization system was converted to 
the software of the higher-precision, large, two-axis 
gimbaled-platform stabilization loop described earlier. The 
only modification required was in the Main Program and the 
Rate-Loop Subprogram for the latter system. The modular 
software design procedure made possible the use of most of 
the building blocks (subroutines) in the implementation of 
the new controller. : 

In general, the study demonstrated the technical 
feasibility of digital control for a wide-bandwidth, high- 
precision type of system. Due to the limited scope of the 
study, the full power of digital control was not utilized, in 
that the control algorithms were constrained by the design 
to emulate their analog prototypes. It is likely that significant 
performance improvements could be achieved by advanced 
control techniques. 

Additional capacity in the TMS32010 remains to 
accommodate improved, more sophisticated compensators. 
Table 3 shows the TMS32010 utilization. 


Table 3. TMS32010 Utilization 
(LOS Stabilization System Rate-Loop) 


Used | Available 
Program 275 4096 
7% 
Memory words words 
Data 76 144 
53% 
Memory words words 
Execution 
Time 


* Based on a 16-MHz (i.e., less than maximum) clock rate. 
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Figure 14. Rate-Loop Normalized Torque-Disturbance Step Response [19] 
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Figure 15. Rate-Loop Normalized Torque-Disturbance Frequency Response [19] 


Other microprocessors that were considered for 
implementing this digital controller system were: Intel 8086, 
Zilog Z-8000, Motorola 68000, and the Fairchild 9445. 
These microprocessors were unable to meet the criterion that 
the maximum allowable time between samples for processing 
be 250 us. Among the signal-processing microprocessors, 
the AMI 2811, while apparently fast enough, has only a 
12 x 12 multiplier; and the Intel 2920 has only four inputs 
and no branching instructions. 

The principal limitation of the TMS32010 was that of 
having eight inputs and eight outputs. Except for this 
restriction, the processor would have been able to carry out 
the processing for both axes of the two-axis gimbaled 
platform. This limitation could be removed by the addition 
of logic circuitry. 
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APPENDIX A 
Development of a Digital Compensator Transfer Function 


The development of a digital equivalent of an analog compensator transfer function using the bilinear transformation with 
frequency prewarping is shown in this appendix. The technique is described in the section, DIGITAL COMPENSATOR 
DESIGN. 


Beginning with an analog prototype transfer function, 


S2 + 68.2 S + 3943 


G(s) = 1000 oy -9512 § + 631x100 


The sampling frequency to be used in converting to a digital equivalent is f = 4020 Hz (i.e., the sampling period 
T, = 1/4020 s = 248.76 x 10-65). 


The characteristic equation of this analog transfer function is: 
S2 + 2512S + 6.31 x 106 = 0 


which fits the standard, second-order form: 


S2 + 2¢unS +2 =0 


The natural frequency w, = \ 6.31x106 = 2511.9713 rad/s 


To compensate for nonlinear mapping of analog-to-digital frequencies by the bilinear transformation method, the natural 
frequency is prewarped according to the formula: 
2 WoT 2 2511.9713 x 248.76 x 10-6 


op Si 2 QAST6x10-o° 2 


= 2597.03 rad/s 
This prewarping scheme matches exactly the natural frequency in the analog and digital domains for the compensator. 


To obtain the prewarped version of the analog transfer function, the complex variable s in the original transfer function is 
replaced with (wo/wp)s. It is therefore convenient to compute the ratio: 


@o _ _2511.9713 _ 
oe OSs 


The prewarped G(s), 1.e., Gp(s) is then computed as: 


1000 0.9672 s)2 + 68.2 (0.9672 s) + 3943 
(0.9672 s)2 + 2512 (0.9672 s) + 6.31 x 106 


Gp(s) 


__s2 + 7051s + 4214.87 __ 
s2 + 2597.16 s+ 6.75 x 106 


1000 


Bilinear transformation is next applied to Gp (s) whereby the continuous variable s is replaced by the expression that involves 
the discrete variable z: 

z—1 

z+1 





. +2 
S = oF 


This produces the discrete transfer function D(z). 


For the compensator, 





D(z) = Gps) a a 
. * T 741 
2 
= 1000 s2_+ 70.51 s+ 4214.87 


s2 + 2597.16 s+ 6.75 x 106 
2 z—1 
248.76x10-6 z+1 





S — 


After further computations, 


- 1.0 — 1.9824 z—1 + 0.9826 z-2 
D@) = 790.76 10 — 1.2548 2-1 + 0.5474 22 

The final step is the gain adjustment in the digital transfer function, This can be accomplished by matching the analog and 
digital gains at some predetermined frequency, for example, DC. 


For the DC case, s = jw = O and from the bilinear transformation: 


. 2st .. \ 
= eee: 





Therefore, at DC, G(O) = D(1). 
For this transfer function, G(O) = 0.6249, D(1) = 0.5072. If G(O) = K x D (1), then the constant K becomes 


0.6249 


0.5072 ~ 1-232! 





The final form.of the digital equivalent transfer function is: 


9824 7-1 + 0.9826 z-2 


— 1041 
: D@) = 810.77 7'9 = 1.2548 2-1 + 0.5474 2-2 


where the gain of 870.77 is the product of K x (the unadjusted digital gain), 1.e., 870.77 = 1.2321 x 706.76. 
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APPENDIX B 
TMS32010 Example Program 


An example TMS32010 program that uses the first- and second-order filter routine to compute several elements of a digital 
control system is provided in this appendix. The program illustrates the concepts of modular software design, data memory 
layout, and cascade implementation of high-order transfer functions. The program was executed on a combination of the 
TMS32010 Emulator (XDS) and Analog Interface Board (AIB) using random noise as input. The input was sampled at a 
4000-Hz rate. 


The following transfer functions are implemented with Q12 scaling: 


0.8352 — 0.2729 z—! + 0.8352 z—2 


900-Hz Notch Filter D(z) = 10 — 0.2729 2-1 + 0.6704 2-2 
| 0.9688 + 1.8341 z—-1 + 0.9688 z—2 
BRO eee ae oy 1.0 + 1.8341 2-1! + 0.9375 2-2 


Other transfer functions (compensators, notch filters) can be implemented in identical fashion by expanding the data structure 
(filter coefficients and states) and making additional filter routine calls to compute these elements. 


The output, as observed on a spectrum analyzer, is shown in Figure B-1. 





Figure B-1. Spectrum Analyzer output (900-Hz and 1800-Hz Notch Filters) 


The first notch from the left is at 900 Hz, the second is at 1800 Hz. The attenuation of the notch frequencies is about 23 
dB in reference to the passband region. 


The program that produced this output is as follows: 
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RIG 


wn 


0001 
0002 
0003 
0004 
O00 
0006 
0007 
0008 
O00? 
0010 
0011 
0012 
0013 
0014 
0015 
0O14 
0017 
0018 
0019 
0070 
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0024 
0027 
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0029 
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0031 
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O035 
0036 
0037 
0036 
OO39 
QO046 
0041 
0042 
0043 
0044 
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0035 
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SY: 


in 


0000 
0000 
0001 


0002 
O003 
0004 
0005 
0006 
0007 
0008 
Q00% 
OOOA 
OOOH 
000C 
0002 
QOOOE 
OOOF 
OO10 


0011 
0012 
0013 
0014 
0015 
0016 
0017 
OO1S 
OO1Y 


OO1A 


OO1B 


OO1IC 


OoOiDd 


OOLE 


QOO000 
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QOOOr 


0000 
0001 


0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
000A 
000B 
o00c 
0o00D 
QOQ0E 
OOOF 
0010 
0010 


0011 
O01 
OOLS 
0014 
OO1S 
OO1L4E 
0017 
0018 
OOLY 
OOL? 


OOLA 
OO1B 


OO1C 
OO1D 
OO1L 
OO1D 
OOLE 


RK ke OO 


SCALE 
* 
# 


MODE 
RATE 
+ 

NOS 
Ni3 
N23 
Dis 
D3 
NOZ 
Ni2 
Nee 
Diz 
a 
NOl 
Nil 
N21 
Dil 
D1 
COEFFS 
# 

XO3 
X13 
X23 
XO2 
X12 
X22 
XO} 
Xil 
X21 
STATES 
+ 
COMAND 
OUTPUT 
# 
MAX14 
MINI4G 
MASK 1 
MASK2 
CINE 

+ 

+ 


DATA 


11°14:45 03-06-85 
FAGE OO01 


~PCS.1 24.107 


IDT “RIG SYS” 
OPTION DUNLST, TUNLST 


The program computes a few sections af the LOS 
Stabilization System Rate Loop. 


Constants 

The constant SCALE relates to the scaling factor 
through a relation: 

scaling factor = 2 ## SCALE (ex. 4094 = 2 ## 12). 
Scale canstants 12 and 15 are available. 

Era 12 SCALING FACTOR (@12) 


Data memory map 
DORG 0 

DATA $ 

DATA $ 
Coefficient table 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


DATA MEM START ADRS 
AIB MODE 
AIB RATE 


FILTER #3 COEFF “S 


FILTER #2 COEFF“S 


DATA 
DATA 
DATA | 
DATA 
DATA 
DATA 
DATA 
EGU 
State va 
DATA 
DATA 
DATA 
LATA 
DATA 
DATA 
DATA 
DATA 
DATA 
Etii 


FILTER #1 COEFF“S 


4A th of Oo oR Ot Ot RR 


i 
pws, 


' table 


~— 


FILTER #: 


het 


STATES 


FILTER #2 STATES 


fl 


FILTER #1 STATES 


tH +h th 4 +R Ot OU Ot OR Oo 


i 
a 


DATA 
DATA 


Sad 


COMMAND INFUT 
SYSTEM OUTPUT 


+ 


DATA $ 
DATA $ 
EG MAX14 
Eau MINI4 
TATA $ 


MAX 2-COMPLEMENT NUM IN 14 BITS 
MIN 2-COMPLEMENT NUM IN 14 BITS 
MASK 

MASK 

ONE 


AUR 0 


RIG SYS 


O0OS8 


OOsF 
0040 
OO61 
0062 
OO0G3 
0064 
0065 
00646 
Q067 
00468 
0069 
0070 
O07 1 
0072 
0073 
0074 
0075 
0076 
0077 


0075 
O07? 
OO8O 


O01 


OOS? 
OOS 
O084 
OO8S 
OOBSG 
0037 
0088 
OOS? 
OO?O 
OO?! 
OOP2 
OO?3 
OOV4 


OOPS 


OOF! 
OOV7 
OOS 
OO? 
OL00 
O1O1 
OL1o02 
OL10O3 
OL104 
OLos 


OLO6G 
OLO7 


0000 
0001 


0004 
OOO 
OO0E 
OOLS 
OO16 
Oo19 
O01C 
OOLE 
0020 


0021 
0022 


0023 
0024 
OO25 
0024 


0027 
0028 
O02? 
002A 


QOO=SR 
O02C 


oo2n 


QOOZE 
QOOZF 
O030 
0031 


OO32 
OO3S 
O034 
OO35 
OOR4 


GO37 
OO3S 
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FOO 
0021 


00072 
QO0O08 
0000 
OnS0 
OF SO 
0000 
0000 
0000 
0000 
7FFF 
0001 
0020 


0021 
FSO0O 
OOR'? 


F400 
OQ27 
FOO 
00223 


0027 
AZ1A 
461A 
7310 
SOLA 


ZLiA 
7010 
7119 
FSOO 
0046 
Foo 
0046 


aC1B 
641B 
71D 
OLR 
4A1B 


F200 


OoORS 


TABLE 


TBLEND 


WAIT 


GET 


B START 

Table in prog memory 
Ei $ 

LATA 2A, >4DB 

DATA 0,0,0,0,0 
DATA 

DATA S968 "7513, 3943, 
DATA 0,0,0 

DATA 0,9,9 

DATA 0.0.0 

DATA o,9 

DATA 32767 ,-32763 
DATA 1 

ERU $—1 


Initialize the system 
Eta $ 
CALL INIT 


Wait on sample 
BIUZ GET 


B WAIT 


Input sample 


Ei $ 

IN COMAND, PAZ 
ZALS COMAND 

XOR MASK 1 

SACL COMAND 
Fracess sample 

LAC COMAND, SCALE 
LARK. ARO, COEFFS 
LARK. AR1, STATES 
CALL FILTR2 
CALL FILTR2 


Dutput sample 


SACH OUTPUT, 14-SCALE 
ZALS NUTPUT 

XUR MASK 2 

SACL QUTRLUT 

CWT MLITFUT | PAZ 


Repeat the sequence 
B WAIT 


POS. 1 34.107 


3421,-1118, 3421,1115,-2744 
-7512, -3340 


esee Yaad Ee, 


11:14:45 O2-04-85 
FAGE OOO? 


RESTART VECTOR 


AIB MODE AND RATE 

FILTER #3 COEFF’S 

900 HZ NOTCH FILTER 
1600 HZ NOTCH FILTER 
FILTER #3 INITIAL STATES 

700 HZ NOTCH INITIAL STATES 

1800 HZ NOTCH INITIAL STATES 
COMMAND INPUT, SYSTEM OUTPUT 

MAX AND MIN 14 BIT NUMBERS 

QINE 


INITIALIZATION ROUTINE 


INPUT COMMAND 

BET COMMAND 
CORRECT A/D FORMAT 
UPDATE COMMAND 


LOAD SCALED COMMAND 


ARO = PTR TO COEFF TABLE 
ARi = FTR TP STATE VAR TABLE 
1200 HZ NOTCH FILTER 


700 HZ NOTCH FILTER 


STORE OUTPUT 

BET OUTPUT 

IORRECT FOR D/A FORMAT 
WPDATE OUTPUT 

AND SEND IT OUT 


60 GET NEXT SAMPLE 


B-3 


RIG SYS 


oO108 
OL10% 
0110 
O111i 
O112 
O113 
0114 
0115 
01146 
0117 
0118 
O119 
0120 
O121 
O122 
0123 
0124 


O125 


0126 


O127 
O125 
O129 
0130 
O131 
O132 
O133 
0134 
O133 
O136 
0137 
0138 
0137 
0140 
O141 
O142 
0143 
0144 
0145 
0144 
0147 
0148 
O149 
0150 
0151 
O1S2 
O133 
O134 
O1355 
0154 
O1357 
0158 
O1S9 
0140 
O161 
0142 
0143 


B-4 


0039 
003A 


0OO3B 
O03 
Oo3D 
OO3E 


OO3F 
0040 
0041 
0042 


0043 
0044 


O045 


0046 


0047 
0048 
004% 
Q04A 
QOO4R 
O04c 


oo04p0 
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0039 
7TFE1 
7FS3B 


6E00 
4850 
JOLIE 
7EZO 
OO3F 
67883 
LOIE 
F400 
OO3F 


4200 
4BO1 


7FS80 


a 
a oO 
Mm 
me PS 


GAz0O 
&éD7 1 
EAD 
4071 
GIP 


bB98 


FAQO 


mt KK x KK OK 


NIT 


XFER 


* 


Th x oR OK KKK 


ILTR2 


xe 


PAGE 0003 


System initialization routine. The routine initializes 
the TMS32010 and other system components. 

The calling sequence is: 

CALL INIT 


Initialize the 32010 


Ell $ 
DINT DISABLE INTERRUPTS 
SOVM SET OVERFLOW MODE 


Initialize Data Memory 


LDPK 0 USE PAGE 0 

LARP 0 USE ARO 

LARK ARO, TBLEND-TABLE INIT PTR TO END OF DATA 
LACK. TBLEND INIT FTR TO END OF TABLE 

EQU $ 

TBLR # XFER FROM PROG TO DATA MEM 

SUB ONE BLIMP PTR DOWN 

BANZ XFER 60 XFER MORE 


Inititialize AIB 


OUT MODE, PAO AIB made 
QUT RATE, PAS AIB RATE 
RET RETURN 


First and second order filter routine. Computes an IIR 
filter using Direct Form II] algorithm and adapts toa a 
scaling scheme defined in the calling program. 


The routine incorporates overflow handling code upon 


computing the intermediate value and the output. 


The calling sequence is: 

ACC = scaled filter input 
ARO = ptr to coeff table 

ARI = ptr to state var table 
CALL FILTR2 


ACC = scaled filter output 

ARO = ptr toa next set of coeff"s 

ARL = ptr to mext set of state vars 
Et $ 

LARP ARI USE AR1L 


Compute intermediate value 


LT #—-, ARO T=X2 

MPY #-, ARI MPY X2#D2 

LTA #+, ARO T=X1, ACC=KU+X2#D2 

MFY #-, ARI MPY X1#D1 

LTA i= T=X2, ACC=KU+X2#D2+X1#D1 
MAR a ARI=FTR TO XO 


Raund, store and check for intermediate overflow 
BLZ LBEL1IO CHECK FOR +/- RESULT 


RIG SY5 


01464 
OL165 
01464 
0167 


OL46 
O16? 
0170 


0171 
O172 
O173 
0174 
O175 


01764 
OL77 
O175 
O179 
0180 
O1sti 
O182 
O153 
01284 
0135 
01386 
0187 
01388 
0159 
0170 
Ovi 


OLS 
O13 
OL194 


O1?S 


OL? 
O197 


O19 
OLY 
O200 
O201 
O202 


0203 
O204 


OFOS 
O2O6 
O207 
OFOS 
O2OD 
O210 
O211 
O712 


OO4E 
OO4F 
OOS50 
0051 
0052 
OO33 
OO34 
OO55 
00546 
0037 


005s 
0034 
OOTA 
OO3B 
OOS 
OosDo 
OOSE 


QOSF 
0060 
0061 
0062 
0063 
0064 
0065 
0064 


0067 
0065 
0067 
QO06A 
QOO4SB 
0040 
0040 
QO6E 
OOGF 
0070 


OO71 
0072 
0073 
0074 
0075 
0074 
0077 
0078 


0079 


Q07A 


NO ERRORS, 
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OoOss 
OBIE 
Sse 
1C1C 
FBOO 
OOSF 
6610 
SOss 
FOO 
OOSF 
0055 
1BiE 


mS 
110 
FOO 
OO3F 
4£41D 


ORS 


LBL10 


+ 
+ 

OOSF LBL2O 

&3AQ 

46D? 1 

7FaY 

B20 

4071 

B70 

4D? 1 

7F SF 


FAOQO 
0071 
OB1E 
SC1B 
1CAic 
FRBOO 
O07? 
SCIC 
FOO 
0078 
0071 
iB1E 
SCLiB 
ici 
FIRDO 
0077 
2CiD 
FOO 
00764 


LBL30O 


+ 
0077 LBL40 
ZLib 
DOTA 
7F SD 


LBL0 


+ 
+ 


NO WARNINGS 


ADD 
SACH 
SUE 
BLEZ 


ZALS 
SACL 
B 


EQU 
SUB 
SACH 
SUB 
BGEZ 


ZALS 
SACL 


Compute 
EQu 

MAR 

MPY 

ZAC 

LTO 

MFY 

LTO 

MFY 
APAC 


ONE, SCALE-1 
#, 16-SCALE 

MAX14, SCALE 
LBL20 


MAX14 
# 


LBLZO 


$ 

ONE, SCALE-1 
#, 16-SCALE 
MINIG, SCALE 
LBL2O 


MINIG 
+ 


filter output 
$ 

#+, ARO 

#- ARI 


#— ARO 
#- ARI 
#-, ARO 
#~ ARI 
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11:14:45 o3-046-85 
FAGE O004 


ROUND 

UPDATE INTERMEDIATE VAL 
SUBTRACT SCALED MAX POS NUMBER 
IF ACC<=0 THEN NO OVERFLOW 


DNVERFLOW, LOAD MAX FOS NUMBER 
UPDATE INTERMEDIATE VALUE 
60, COMPUTE OUTPUT 


ROUND 

WIFDATE INTERMEDIATE VAL 
SUBTRACT SCALED MIN NEG NUMBER 
IF ACC>=0 THEN NO OVERFLOW 


DVERFLOW, LOAD MIN NEG NUMBER 
UPDATE INTERMEDIATE VALUE 


WISE ARQ 


MPY X2Z#N2 
CLR ACC 
T=X1, ACC=X2eNZ, UPDATE X2 


MEY Xi#N1 

T=XO, ACC=X2eNZ+X1#N1, UPDATE X1 
MEY XO#NO 

ACC=X2#N2+X 14#N1+XO#NO 


Check for output averflow 


BL.Z 


ADD 
SALH 
SUB 
BLEZ 


LAC 
B 


Ei 
SUB 
SACH 
SUB 
BGEZ 


LAC 
B 


EG 
LAL 
Ett 
RET 


END 


LBL:30 


ONE, SCALE-1 


DUTPUT , 14-SCALE 


MAX14, SCALE 
LBL4o 


MAX14, SCALE 
LBLSO 


$ 
ONE, SCALE-1 


OUTPUT, 14-SCALE 


MINi4, SCALE 
LBL40 


MINI4, SCALE 
LBLSO 


$ 
OUTPUT, SCALE 
$ 


CHECK FOR +/- RESULT 


ROUND 

WIPDATE OUTPUT 

SUBTRACT SCALED MAX FOS NUMBER 
IF ACC<=0 THEN NO OVERFLOW 


OVERFLOW, LOAD MAX FOS NUMBER 
bO. RETURN 


wie gs 


ROLIND 

LIFDATE OUTPUT 

SUBTRACT SCALED MIN NEG NUMBER 
IF ACC>=0 THEN NO OVERFLOW 


OVERFLOW, LOAD MIN NEG NUMBER 
GO, RETURN 


ent gy 


RESTORE ACC 


RETURN 


TI Worldwide 
Sales Offices 


ALABAMA: Huntsville: 500 Wynn Drive, Suite 514, 
Huntsville, AL 35805, (205) 837-7530. 


ARIZONA: Phoenix: 8825 N. 23rd Ave., Phoenix, 
AZ 85021, (602) 995-1007. 


CALIFORNIA: Irvine: 17891 Cartwright Rd., Irvine, 
CA 92714, (714) 660-1200; Sacramento: 1900 Point 
West Way, Suite 171, Sacramento, CA 95815, 

(916) 929-1521; San Diego: 4333 View eke Ave., 
Suite B., San Diego, CA 92123, (619) 278-9601: 
Santa Clara: 5353 Betsy Ross Dr., Santa Clara, CA 
95054, (408) 980-9000; Torrance: 19505 Hamilton St., 
Bidg. A, Suite 1, Torrance, CA 90502, (213) 217-7010; 
Woodland Hills: 21220 Erwin St., Woodland Hills, 
CA 91367, (213) 704-7759. 


COLORADO: Aurora: 1400 S. Potomac Ave., 
Suite 101, Aurora, CO 80012, (303) 368-8000. 


CONNECTICUT: Wallingford: 9 Barnes Industrial 
Park Rd., Barnes Industrial Park, Wallingford, 
CYT 06492, (203) 269-0074. 


FLORIDA: Ft. Lauderdale: 2765 N.W. 62nd St., 

Ft. Lauderdale, FL 33309, (305) 973-8502; Maitland: 
2601 Maitland Center Parkway, Maitland, FL 32751, 
(305) 660-4600; Tampa: 5010 W. Kennedy Bilvd., 
Suite 101, Tampa, FL 33609, (813) 870-6420. 


GEORGIA: Norcross: 5515 Spalding Drive, Norcross, 
GA 30092, (404) 662-7900 


ILLINOIS: Arlington Heights: 515 W. Algonquin, 
Arlington Heights, IL 60005, (312) 640-2925. 


INDIANA: Ft. Wayne: 2020 Inwood Dr., Ft. Wayne, 
IN 46815, (219) 424-5174; Indianapolis: 2346 S. 
Lynhurst, Suite J-400, Indianapolis, IN 46241, 
(317) 248-8555. 


1OWA: Cedar Rapids: 373 Collins Rd. NE, Suite 200, 
Cedar Rapids, !A 52402, (319) 395-9550. 


MARYLAND: Baltimore: 1 Rutherford Pl., 
7133 Rutherford Rd., Baltimore, MD 21207, 
(301) 944-8600. 


MASSACHUSETTS: Waitham: 504 Totten Pond Rad., 
Waltham, MA 02154, (617) 895-9100. 


MICHIGAN: Farmington Hills: 33737 W. 12 Mile Rd., 
Farmington Hills, M! 48018, (313) 553-1500. 


MINNESOTA: Eden Prairie: 11000 W. 78th St., 
Eden Prairie, MN 55344 (612) 828-9300. 


MISSOURI: Kansas City: 8080 Ward Pkwy., Kansas 
City, MO 64114, (816) 523-2500; St. Louis: 

11861 Westline Industrial Drive, St. Louis, 

MO 63141, (314) 569-7600. 


NEW JERSEY: Iselin: 485E U.S. Route 1 South, 
Parkway Towers, Iselin, NJ 08830 (201) 750-1050 


NEW MEXICO: Albuquerque: 2820-D Broadbent Pkwy 
NE, Albuquerque, NM 87107, (505) 345-2555. 


NEW YORK: East Syracuse: 6365 Collamer Dr., East 
Syracuse, NY 13057, (315) 463-9291; Endicott: 112 
Nanticoke Ave., P.O. Box 618, Endicott, 

NY 13760, (607) 754-3900; Melville: 1 Huntington 
Quadrangle, Suite 3C10, P.O. Box 2936, Melville, 
NY 11747, (516) 454-6600; Pittsford: 2851 Clover St., 
Pittsford, NY 14534, (716) 385-6770; Poughkeepsie: 
385 South Rd., Poughkeepsie, NY 12601, 

(914) 473-2900. 


NORTH CAROLINA: Charlotte: 8 Woodlawn Green, 
Woodlawn Rd., Charlotte, NC 28210, (704) 527-0930; 
Raleigh: 2809 Highwoods Bivd., Suite 100, Raleigh, 
NC 27625, (919) 876-2725. 


OHIO: Beachwood: 23408 Commerce Park Rad., 
Beachwood, OH 44122, (216) 464-6100; Dayton: 
Kingsley Bldg., 4124 Linden Ave., Dayton, OH 45432, 
(513) 258-3877. 


OKLAHOMA: Tulsa: 7615 East 63rd Place, 
3 Memorial Place, Tulsa, OK 74133, (918) 250-0633. 


OREGON: Beaverton: 6700 SW 105th St., Suite 110, 
Beaverton, OR 97005, (503) 643-6758. 


PENNSYLVANIA: Ft. Washington: 260 New York Dr., 
Ft. Washington, PA 19034, (215) 643-6450; 
Coraopolis: 420 Rouser Rd., 3 Airport Office Park, 
Coraopolis, PA 15108, (412) 771-8550. 


PUERTO RICO: Hato Rey: Mercantil Plaza Bidg., 
Suite 505, Hato Rey, PR 00919, (809) 753-8700. 


TEXAS: Austin: 12501 Research Bivd., 

P.O. Box 2909, Austin, TX 78723, (512) 250-7655; 
Richardson: 1001 E. Campbell Rd., 

Richardson, TX 75080, 

(214) 680-5082; Houston: 9100 Southwest Frwy., 
Suite 237, Houston, TX 77036, (713) 778-6592: 
San Antonio: 1000 Central Parkway South, 

San Antonio, TX 78232, (512) 496-1779. 


UTAH: Murray: 5201 South Green SE, Suite 200, 
Murray, UT 84107, (801) 266-8972. 


VIRGINIA: Fairfax: 3001 Prosperity, Fairfax, VA 
22031, (703) 849-1400. 


WASHINGTON: Redmond: 5010 148th NE, Bldg B, 
Suite 107, Redmond, WA 98052, (206) 881-3080. 


WISCONSIN: Brookfield: 450 N. Sunny Slope, 
Suite 150, Brookfield, WI 53005, (414) 785-7140. 


CANADA: Nepean: 301 Moodie Drive, Mallorn 
Center, Nepean, Ontario, Canada, K2H9C4, 

(613) 726-1970. Richmond Hill: 280 Centre St. E., 
Richmond Hill L4C1B1, Ontario, Canada 

(416) 884-9181; St. Laurent: Ville St. Laurent Quebec, 
9460 Trans Canada Hwy., St. Laurent, Quebec, 
Canada H4S1R7, (514) 334-3635. 


ARGENTINA: Texas Instruments Argentina 
S.A.1.C.F.: Esmeraida 130, 15th Floor, 1035 Buenos 
Aires, Argentina, 1+ 394-3008. 


AUSTRALIA (& NEW ZEALAND): Texas Instruments 
Australia Ltd.: 6-10 Talavera Rd., North Ryde 
(Sydney), New South Wales, Australia 2113, 

2 + 887-1122; 5th Floor, 418 St. Kilda Road, 
Melbourne, Victoria, Australia 3004, 3 + 267-4677; 
171 Philip Highway, Elizabeth, South Australia 5112, 
8 + 255-2066. 


AUSTRIA: Texas Instruments Ges.m.b.H.: 
Industriestrabe B/16, A-2345 Brunn/Gebirge, 
2236-846210. 


BELGIUM: Texas Instruments N.V. Belgium S.A.: 
Mercure Centre, Raketstraat 100, Rue de la Fusee, 
1130 Brussels, Belgium, 2/720.80.00. 


BRAZIL: Texas Instruments Electronicos do Brasil 
Ltda.: Rua Paes Leme, 524-7 Andar Pinheiros, 05424 
Sao Paulo, Brazil, 0815-6166. 


DENMARK: Texas Instruments A/S, Mairelundvej 
46E, DK-2730 Herlev, Denmark, 2 - 91 74 00. 


FINLAND: Texas Instruments Finland OY: 
Teollisuuskatu 19D 00511 Helsinki 51, Finland, (90) 
701-3133. 


FRANCE: Texas Instruments France: Headquarters 
and Prod. Plant, BP 05, 06270 Villeneuve-Loubet, 
(93) 20-01-01; Paris Office, BP 67 8-10 Avenue 
Morane-Sauinier, 78141 Velizy-Villacoublay, 

(3) 946-97-12; Lyon Sales Office, L’Oree D’Ecully, 
Batiment B, Chemin de la Forestiere, 69130 Ecully, 
(7) 833-04-40; Strasbourg Sales Office, Le Sebastopol 
3, Quai Kleber, 67055 Strasbourg Cedex, 

(88) 22-12-66; Rennes, 23-25 Rue du Puits Mauger, 
35100 Rennes, (99) 31-54-86; Toulouse Sales Office, 


Le Peripole—2, Chemin du Pigeonnier de la Cepiere, 


31100 Toulouse, (61) 44-18-19; Marseille Sales Office, 
Noilly Paradis—146 Rue Paradis, 13006 Marseille, 
(91) 37-25-30. 


TEXAS 
INSTRUMENTS 


Creating useful products 
and services for you 


GERMANY (Fed. Republic of Germany): Texas 
Instruments Deutschland GmbH: Haggertystrasse 1, 
D-8050 Freising, 8161 + 80-4591; Kurfuerstendamm 
195/196, D-1000 Berlin 15, 30 + 882-7365; Ill, Hagen 
43/Kibbelstrasse, .19, D-4300 Essen, 201-24250; 
Frankfurter Allee 6-8, D-6236 Eschborm 1, 
06196 + 8070; Hamburgerstrasse 11, D-2000 Hamburg 
76, 040 + 220-1154, Kirchhorsterstrasse 2, D-3000 
Hannover 51, 511+648021; Maybachstrabe 11, 
D-7302 Ostfildern 2-Nelingen, 711 +547001; 
Mixikoring 19, D-2000 Hamburg 60, 40 + 637 +0061; 
Postfach 1309, Roonstrasse 16, D-5400 Koblenz, 

261 + 35044. 


HONG KONG (+ PEOPLES REPUBLIC OF CHINA): 
Texas Instruments Asia Ltd., 8th Floor, World 
Shipping Ctr., Harbour City, 7 Canton Rd., Kowloon, 
Hong Kong, 3 + 722-1223. 


IRELAND: Texas Instruments (Ireland) Limited: 
Tees Rd., Stillorgan, County Dublin, Eire, 
1 831311. 


ITALY: Texas Instruments Semiconduttori Italia Spa: 
Viale Delle Scienze, 1, 02015 Cittaducale (Rieti), 
Italy, 746 694.1; Via Salaria KM 24 (Palazzo Cosma), 
Monterotondo Scalo (Rome), Italy, 6 + 9003241; Viale 
Europa, 38-44, 20093 Cologno Monzese (Milano), 

2 2532541; Corso Svizzera, 185, 10100 Torino, Italy, 
aaa Via J. Barozzi 6, 40100 Bologna, Italy, 51 

1. 


JAPAN: Texas Instruments Asia Ltd.: 4F Aoyama 
Fuji Bldg., 6-12, Kita Aoyama 3-Chome, Minato-ku, 
Tokyo, Japan 107, 3-498-2111; Osaka Branch, 5F, 
Nissho Iwai Bldg., 30 Imabashi 3- Chome, 
Higashi-ku, Osaka, Japan 541, 06-204-1881; Nagoya 
Branch, 7F Daini Toyota West Blidg., 10-27, Meieki 
4-Chome, Nakamura-ku Nagoya, Japan 

450, 52-583-8691. 


KOREA: Texas Instruments Supply Co.: 3rd Floor, 
Samon Bidg., Yuksam-Dong, Gangnam-ku, 
135 Seoul, Korea, 2 + 462-8001. 


MEXICO: Texas Instruments de Mexico S.A.: Mexico 
City, AV Reforma No. 450 — 10th Floor, Mexico, 
D.F., 06600, 5+ 514-3003. 


MIDDLE EAST: Texas Instruments: No. 13, 1st Floor 
Mannai Bldg., Diplomatic Area, P.O. Box 26335, 
Manama Bahrain, Arabian Gulf, 973 + 274681. 


NETHERLANDS: Texas Instruments Holland B.V., 
P.O. Box 12995, (Bullewijk) 1100 CB Amsterdam, 
Zuid-Oost, Holland 20 + 5602911. 


NORWAY: Texas Instruments Norway A/S: PB106, 
Refstad 131, Osio 1, Norway, (2) 155090. 


PHILIPPINES: Texas Instruments Asia Ltd.: 14th 
Floor, Ba- Lepanto Bidg., 8747 Paseo de Roxas, 
Makati, Metro Manila, Philippines, 2 + 8188987. 


PORTUGAL: Texas Instruments Equipamento 
Electronico (Portugal), Lda.: Rua Eng. Frederico 
Ulrich, 2650 Moreira Da Maia, 4470 Maia, Portugal, 
2-948-1003. 


SINGAPORE (+ INDIA, INDONESIA, MALAYSIA, 
THAILAND): Texas Instruments Asia Ltd.: 12 Lorong 
Bakar Batu, Unit 01-02, Kolam Ayer Industrial Estate, 
Republic of Singapore, 747-2255. 


SPAIN: Texas Instruments Espana, S.A.: C/Jose 
Lazaro Galdiano No. 6, Madrid 16, 1/458.14.58. 


SWEDEN: Texas Instruments International Trade 
Corporation (Sverigefilialen): Box 39103, 10054 
Stockholm, Sweden, 8 - 235480. 


SWITZERLAND: Texas Instruments, Inc., Reidstrasse 
6, CH-8953 Dietikon (Zuerich) Switzerland, 
1-740 2220. 


TAIWAN: Texas Instruments Supply Co.: Room 903, 
205 Tun Hwan Rd., 71 Sung-Kiang Road, Taipei, 
Taiwan, Republic of China, 2 + 521-9321. 


UNITED KINGDOM: Texas Instruments Limited: 
Manton Lane, Bedford, MK41 7PA, England, 0234 
67466; St. James House, Wellington Road North, 
Stockport, SK4 2RT, England, 61+ 442-7162. BK 





TI Sales Offices 


ALABAMA: Huntsville (205) 837-7530. 
ARIZONA: Phoenix (602) 995-1007. 


CALIFORNIA: Irvine (714) 660-1200; 
Sacramento (916) 929-1521; 

San Diego (619) 278-9601; 

Santa Clara (408) 980-9000; 
Torrance (213) 217-7010; 

Woodland Hills (213) 704-7759. 


COLORADO: Aurora (303) 368-8000. 
CONNECTICUT: Wallingford (203) 269-0074. 


FLORIDA: Ft. Lauderdale (305) 973-8502; 
Maitland (305) 660-4600; Tampa (813) 870-6420. 


GEORGIA: Norcross (404) 662-7900. 
ILLINOIS: Arlington Heights (312) 640-2925. 


INDIANA: Ft. Wayne (219) 424-5174; 
Indianapolis (317) 248-8555. 


IOWA: Cedar Rapids (319) 395-9550. 
MARYLAND: Baltimore (301) 944-8600. 
MASSACHUSETTS: Waltham (617) 895-9100. 
‘MICHIGAN: Farmington Hills (313) 553-1500. 
MINNESOTA: Eden Prairie (612) 828-9300. 


MISSOURI: Kansas City (816) 523-2500; 
St. Louis (314) 569-7600. 


NEW JERSEY: Iselin (201) 750-1050. 
NEW MEXICO: Albuquerque (505) 345-2555. 


NEW YORK: East Syracuse (315) 463-9291; 
Endicott (607) 754-3900; Melville (516) 454-6600; 
Pittsford (716) 385-6770; 

Poughkeepsie (914) 473-2900. 


NORTH CAROLINA: Charlotte (704) 527-0930; 
Raleigh (919) 876-2725. 


OHIO: Beachwood (216) 464-6100; 
Dayton (513) 258-3877. 


OKLAHOMA: Tulsa (918) 250-0633. 
OREGON: Beaverton (503) 643-6758. 


PENNSYLVANIA: Ft. Washington (215) 643-6450; 
Coraopolis (412) 771-8550. 


PUERTO RICO: Hato Rey (809) 753-8700 


TEXAS: Austin (512) 250-7655; 
Houston (713) 778-6592; Richardson (214) 680-5082; 
San Antonio (512) 496-1779. 


UTAH: Murray (801) 266-8972. 

VIRGINIA: Fairfax (703) 849-1400. 
WASHINGTON: Redmond (206) 881-3080. 
WISCONSIN: Brookfield (414) 785-7140. 


CANADA: Nepean, Ontario (613) 726-1970; 
Richmond Hill, Ontario (416) 884-9181; 
St. Laurent, Quebec (514) 334-3635. 


TI Regional 
Technology Centers 


CALIFORNIA: Irvine (714) 660-8140, 
Santa Clara (408) 748-2220. 


GEORGIA: Norcross (404) 662-7945. 
ILLINOIS: Arlington Heights (312) 640-2909. 
MASSACHUSETTS: Waltham (617) 890-6671. 
TEXAS: Richardson (214) 680-5066. 


TI Technical 
Support Center 


TOLL FREE: (800) 232-3200, U.S., 
(214) 995-2171, Canada 


TI Distributors 


Cae DISTRIBUTORS IN 
USA 

Arrow Electronics 

Diplomat Electronics 

ESCO Electronics 

General Radio Supply Company 

Graham Electronics 

Harrison Equipment Co. 

international Electronics 

JACO Electronics 

Kierulff Electronics 

LCOMP, Incorporated 

Marshall Industries 

Milgray Electronics 

Newark Electronics 

Rochester Radio Supply 

Time Electronics 

R.V. Weatherford Co. 

Wyle Laboratories 


Ti AUTHORIZED DISTRIBUTORS IN 
CANADA 

CESCO Electronics, Inc. 

Future Electronics 

ITT Components 

L.A. Varah, Ltd. 


ALABAMA: Arrow (205) 882-2730; 
Kierulff (205) 883-6070; Marshall (205) 881-9235. 


ARIZONA: Arrow (602) 968-4800; 
Kierulff (602) 243-4101; Marshall (602) 968-6181; 
Wyle (602) 249-2232. 


CALIFORNIA: Los Angeles/Orange County: 

Arrow (818) 701-7500, (714) 838-5422; 

Kierulff (213) 725- 0325, (714) 731 S71 (714) 220-6300; 
Marshall (213) 999- 5001, (818) 442- 204, 

(714) 660-0951; R.V. Weatherford via 634-9600, 
(213) 849-3451, (714) 623-1261; Wyle (213) 322-8100, 
(714) 863-9953; Sacramento: Arrow (916) 925-7456; 
Wyle (916) 638-5282; San Diego: 

Arrow (619) 565-4800; Kierulff (619) 278-2112; 
Marshall (619) 578-9600; Wyle (619) 565-9171: 

San Francisco Bay Area: Arrow (408) 745-6600; 
(415) 487-4600; Kierulff (408) 971-2600; 

Marshall (408) 732-1100; Wyle (408) 727-2500; 
Santa Barbara: R.V. Weatherford (805) 965-8551. 


COLORADO: Arrow (303) 696-1111; 
Kierulff (303) 790-4444; Wyle (303) 457-9953. 


CONNECTICUT: Arrow (203) 265-7741; 
Diplomat (203) 797-9674; Kierulff (203) 265-1115; 
Marshall (203) 265-3822; Milgray (203) 795-0714. 


FLORIDA: Ft. Lauderdale: Arrow (305) 429-8200; 
Diplomat (305) 974-8700; ee (305) 486-4004; 
Orlando: Arrow (305) 725-1480; 

Milgray (305) 647-5747; Tampe. 

Arrow (813) 576-8995; Diplomat (813) 443-4514; 
Kierulff (813) 576-1966. 


GEORGIA: Arrow (404) 449-8252; 
Kierulff (404) 447-5252; Marshall (404) 923- 5750. 
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ILLINOIS: Arrow (312) 397-3440; 
Diplomat (312) 595-1000; Kierulff (312) 640-0200; 
Marshall (312) 490-0155; Newark (312) 638-4411. 


INDIANA: Indianapolis: Arrow (317) 243-9353; 
Graham (317) 634-8202; 
Ft. Wayne: Graham (219) 423-3422. 


IOWA: Arrow (319) 395-7230. 


KANSAS: Kansas City: Marshall (913) 492-3121; 
Wichita: LCOMP (316) 265-9507. 


MARYLAND: Arrow (301) 247-5200; 
Diplomat (301) 995-1226; Kierul ff (301) 636-5800; 
Milgray (301) 793-3993. 


MASSACHUSETTS: Arrow (617) 933-8130; 
Diplomat (617) 935-6611; Kierulff (617) 667-8331; 
Marshall (617) 272-8200; Time (617) 935-8080. 


MICHIGAN: Detroit: Arrow (313) 971-8220; 
Marshall (313) 525-5850; Newark (313) 967-0600; 
Grand Rapids: Arrow (616) 243-0912. 


MINNESOTA: Arrow (612) 830-1800; 
Kierulff (612) 941-7500; Marshall (612) 559-2211. 


MISSOURI: Kansas City: LCOMP (816) 221-2400; 
St. Louis: Arrow (314) 567-6888; 
Kierulff (314) 739-0855. 


NEW HAMPSHIRE: Arrow (603) 668-6968. 


NEW JERSEY: Arrow (201) 575-5300, (609) 596-8000; 
Diplomat (201) 785-1830; 

General Radio (609) 964-8560; Kierulff (201) 575-6750; 
(609) 235-1444; Marshall (201) 882-0320, 

(609) 234-9100; Milgray (609) 983-5010. 


NEW MEXICO: Arrow (505) 243-4566; 
International Electronics (505) 345-8127. 


NEW YORK: Long Island: Arrow (516) 231-1000; 
Diplomat (516) 454-6400; JACO (516) 273-5500; 
Marshall (516) 273-2053; ey (516) 420- 9800; 
Rochester: Arrow (716) '275- 0300; 

Marshall (716) 235-7620; 

Rochester Radio Supply (716) 454-7800; 
Syracuse: Arrow (315) 652-1000; 

Diplomat (315) 652-5000; Marshall (607) 754-1570. 


NORTH CAROLINA: Arrow (919) 876-3132, 
(919) 725-8711; Kierulff (919) 872-8410. 


OHIO: Cincinnati: Graham (513) 772-1661; 
Cleveland: Arrow (216) 248-3990; 

Kierulff (216) 587-6558; Marshall (216) 248-1788. 
Columbus: Graham (614) 895-1590; 

Dayton: Arrow (513) 435-5563; 

ESCO (513) 226-1133; Kierulff (513) 439-0045; 
Marshall (513) 236- 8088. 


OKLAHOMA: Arrow (918) 665-7700; 
Kierulff (918) 252-7537. 


OREGON: Arrow (503) 684-1690; 
Wyle (503) 640-6000; Marshall (503) 644-5050. 


PENNSYLVANIA: Arrow (412) 856-7000, 
(215) 928-1800; General Radio (215) 922-7037. 


TEXAS: Austin: Arrow (512) 835-4180; 

Kierulff (612) 835-2090; Marshall (512) 458-5654; 
Wyle (512) 834-9957; Dallas: Arrow (214) 380-6464; 
International Electronics (214) 233-9323; 

Kierulff (214) 343-2400; Marshall (214) 233-5200; 
Wyle (214) 235-9953; 

El Paso: International Electronics (915) 598-3406; 
Houston: Arrow (713) 530-4700; 

Marshall (713) 789-6600; 

Harrison Equipment (713) 879-2600; 

Kierulff (713) 530-7030; Wyle (713) 879-9953. 


UTAH: Diplomat (801) 486-4134; 
Kierulff (801) 973-6913; Wyle (801) 974-9953. 


VIRGINIA: Arrow (804) 282-0413. 


WASHINGTON: Arrow (206) 643-4800; 
Kierulff (206) 575-4420; Wyle (206) 453-8300; Marshall 
(206) 747-9100. 


WISCONSIN: Arrow (414) 764-6600; 
Kierulff (414) 784-8160. 


CANADA: Calgary: Future (403) 259-6408; Varah 
(403) 230-1235; Edmonton: Future (403) 486-0974; 
Montreal: CESCO (514) 735-5511; Future 

(514) 694-7710; ITT Components (514) 735-1177; 
Ottawa: CESCO (613) 226-6903; Future 

(613) 820-8313; ITT Components (613) 226-7406; 
Varah (613) 726-8884; Quebec City: CESCO 
(418) 687-4231; Toronto: CESCO (416) 661-0220; 
Future (416) 638-4771; ITT Components 

(416) 630-7971; Varah (416) 516-9311; 
Vancouver: Future (604) 438-5545; Varah 

(604) 873-3211; ITT Components (604) 270-7805; 
Winnipeg: Varah (204) 633-6190 
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